Lighting Condition Analysis for Mars’ Moon Phobos 
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Abstract—This study used high fidelity computer simulation to 
investigate the lighting conditions, specifically the solar radia- 
tion flux over the surface, on Phobos. Ephemeris data from 
the Jet Propulsion Laboratory (JPL) DE405 model was used to 
model the state of the Sun, Earth, Moon, and Mars. An occul- 
tation model was developed to simulate Phobos’ self-shadowing 
and its solar eclipses by Mars. The propagated Phobos state was 
compared with data from JPL’s Horizon system to ensure the ac- 
curacy of the result. Results for Phobos lighting conditions over 
one Martian year are presented, which include the duration of 
solar eclipses, average solar radiation intensity, surface exposure 
time, available energy per unit area for sun tracking arrays, and 
available energy per unit area for fixed arrays (constrained by 
incident angle). The results show that: Phobos’ solar eclipse 
time varies throughout the Martian year, with longer eclipse 
durations during the Martian spring and fall seasons and no 
eclipses during the Martian summer and winter seasons; solar 
radiation intensity is close to minimum at the summer solstice 
and close to maximum at the winter solstice; exposure time per 
orbit is relatively constant over the surface during the spring 
and fall but varies with latitude during the summer and winter; 
and Sun tracking solar arrays generate more energy than a fixed 
solar array. A usage example of the result is also present in this 
paper to demonstrate the utility. 
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1. INTRODUCTION 


Many studies have suggested that a human mission to Phobos 
is an important precursor and catalyst for NASA’s human 
missions to Mars [1], [2], [3], [4]. Since Phobos does not 
have an atmosphere, the mission does not require a complex 
aerocapture system for surface exploration. Moreover, a 
mission to Phobos can be achieved with smaller minimum 
AV than a mission to Mars [2], [3]. Many technologies 
developed for a human Phobos mission are fundamental 
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elements for a successful human Mars mission[3], [4]. Addi- 
tionally, tele-operating Mars surface robots from Phobos will 
incur significantly lower communication delays compared to 
tele-operating from Earth. This near real time commanding 
could allow an in-depth study of Mars and potentially allow 
robotic constructions to prepare for a later human arrival 
[2], [4]. Besides the technical benefits, a human mission to 
Phobos will stimulate public interest in space exploration, 
inspiring the next generation of scientists, engineers, and 
space explorers [2], [3], [4]. 


To ensure the success of a Phobos mission, a comprehensive 
understanding of Phobos’ environment is required. Solar ra- 
diation plays a crucial role in power and thermal subsystems 
for a solar-powered vehicle. A thorough understanding of so- 
lar radiation on Phobos allows engineers to appropriately size 
the solar arrays and the batteries for the power subsystem. It 
also provides important information for astronauts to choose 
the optimal path of the mobility vehicle. This paper focuses 
on the power subsystem as an example to demonstrate the 
usage of Phobos lighting condition data. 


This study uses high fidelity computer simulation to inves- 
tigate the lighting conditions, specifically solar radiation, on 
the Phobos surface over one Martian year. The computer sim- 
ulation was developed using Johnson Space Center’s (JSC) 
in-house simulation tools in order to: (a) model the states of 
the Sun, Earth, Mars, and the Moon using JPL DE405 model; 
(b) model the orbit of Phobos, its surface, and its gravitational 
field; (c) model the occultation of Phobos’ surface due to 
solar eclipse by Mars and self-shadowing. 


The result of the analysis contains average radiation intensity 
on Phobos and surface radiation flux. This data allows us 
to determine surface exposure time based on Phobos orbit, 
duration of solar eclipse by Mars, radiation flux for a fixed 
solar array, and radiation flux for a perfectly sun-tracking 
solar array. The data showed that the lighting condition 
on Phobos’ surface changes dramatically from one Martian 
season to another. This presents a challenge in the vehicle 
design for a long duration mission on Phobos. 


The remainder of this paper is organized as follows. Sec- 
tion 2 describes development of the high fidelity computer 
simulation. Section 3 presents the lighting condition analysis 
result of Phobos for one Martian year. Section 4 shows the 
utility of the lighting condition data. Section 5 enumerates 
the conclusions of this study. 


2. SIMULATION DEVELOPMENT 


A high fidelity model was developed using JSC in-house 
simulation tools to generate the data needed for the analysis. 
The simulation contains three major components including 
an environment model, a Phobos model, and an occultation 
model. These three components are inter-connected to pro- 
duce the presented results. Figure | illustrates an overview of 
the simulation structure. 
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Figure 1. Phobos lighting simulation structure. 


Environment Model 


The environment model, developed using JSC Engineering 
Orbital Dynamic (JEOD), contains models of the Sun, Earth, 
Moon, and Mars. The state data of these planets and the 
Sun were generated using JPL ephemeris file DE405, which 
contains highly accurate ephemeris data of major bodies in 
our solar system[5]. A spherical gravitational model was 
used for the Sun, Earth, and Moon. A 10-degree spherical 
harmonic gravitational model was used for Mars. JEOD 
manages the state of all bodies, including Phobos, using 
a tree-structured reference frame system that calculates an 
accurate gravitational acceleration for Phobos’ state propa- 
gation. Deimos was not included in the simulation due to its 
negligible gravitational effect. 


Phobos Model 


The Phobos model contains a dynamics model, a surface 
model, and a gravitational model. The dynamics model 
was managed by JEOD for state propagation to obtain an 
accurate state of Phobos with respect to other bodies. The 
surface model simulates the detailed shape of Phobos, which 
was used in the gravitational and the occultation models. 
A polyhedron gravitational representation technique is used 
to determine the gravitational field of Phobos. Figure 2 
illustrates one of the STL models of Phobos used in the study. 


Occultation Model 


This study used an occultation model that determines if any 
STL facet is being shadowed by any other facet representing 
the surface of Phobos. The occultation model uses a point 
to point ray-trace from the subject facet to a point source, in 
this case the Sun. The model then solves for any intersection 
of this ray-trace by all other facets, which results in the 
subject facet being shadowed if it is intersected. An optimiza- 
tion routine groups multiple facets into spherical coordinate 
height fields so that only facets of relevant height fields are 
considered during the occultation determination. 





Figure 2. Phobos STL model used in the study. 


3. ANALYSIS RESULTS 


The lighting condition of Phobos’ surface from July 05, 
2030 to May 22, 2032 was generated using the simulation 
developed for the study. The time span covers all four Martian 
seasons, and Table | shows equinoxes and solstices of Mars 
using the northern hemisphere as a reference[6]. 


The propagated state of Phobos was validated with data from 
JPL Horizon system, and Figure 3 shows the position and 
velocity error of the propagated Phobos state. The average 
position error is 26.7m with a standard deviation of 21.6m. 
The average velocity error is 0.006044m/s with a standard 
deviation of 0.0044m/s. 


Figure 4 shows the instantaneous lighting condition of Pho- 
bos on April 01, 2030 00:00:00 using a 5040 facet STL 
model. This date is between a Martian winter and spring 
such that the solar flux is higher in the southern hemisphere. 
The geolocical map of Phobos used in Figure 4, and later 
in Figure 6 through Figure 8 was obtained from the U.S. 
Geological Survey (USGS)[7]. 


The following subsections present post-processed results us- 
ing the lighting conditions data. The equations used for 
computing the results are presented in the Appendix. 


Average Radiation Intensity and Solar Eclipse 


Figure 5 shows the average radiation intensity on Phobos and 
solar eclipse time per orbit over one Martian year. 


The variation in solar radiation intensity on Phobos is due 
to its distance from the Sun. The minimum intensity occurs 
during the summer as Mars is at its maximum distance from 
the Sun, whereas the maximum intensity occurs during the 
Martian winter as Mars is at its minimum distance from the 
Sun. Due to the high inclination angle of Phobos’ orbit with 
respect to Mars’ orbital plane (about 24.19°), the eclipse time 
ranges from 55 minutes at the spring and fall equinoxes to 
zero minutes at the summer and winter solstices. 


Surface Exposure Time per Orbit 


Using the lighting conditions data, the sunlight exposure time 
of Phobos’ surface per orbital period was determined. The 


Table 1. Martian season (northern hemisphere perspective) 
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Figure 3. Propagated state error of Phobos 


results show that: (a) exposure time 1s relatively uniform over 
the surface at the spring equinox with a shorter exposure at 
the sub-Mars area (surface facing Mars) and a longer expo- 
sure at the anti-Mars area; (b) moving towards the summer 
solstice, exposure decreases in the southern hemisphere and 
increases in the northern hemisphere, reaching a maximum 
variance at the summer solstice; (c) exposure time becomes 
relatively uniform again as time progresses from summer 
to fall; (d) then finally, exposure increases in the southern 
hemisphere and decreases in the northern hemisphere as time 
progresses towards the winter solstice. Figure 6 shows four 
instances of the surface exposure time per Phobos orbit at the 
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Figure 4. Lighting condition on April 01, 2030 00:00:00 
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Figure 5. Radiation intensity and solar eclipse of Phobos 


equinoxes and solstices. This information may prove useful 
in determining the sizing requirement of the power subsystem 
of a lander such that it may operate for extended periods, 
including periods of shadowing. 


Radiation Energy per Orbit 


The radiation energy results contain two parts: (a) the avail- 
able energy per unit area for a fixed solar array; (b) the 
available energy per unit area for a perfectly sun tracking 
array. Based on these results, a sun tracking solar array could 
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Figure 6. Surface exposure time per Phobos’ orbit at equinoxes and solstices 


generate 50 to 75 percent more energy then a fixed solar array, 
depending on the time of the Martian year and the location on 
Phobos. Both cases measure the amount of solar energy per 
unit area that is available for one Phobos orbit as a function 
of time (100% efficient solar array). 


Available energy for Fixed Solar Array—tepresents a conser- 
vative scenario for energy generation. It assumed the solar 
array is orientated parallel to the surface of Phobos without 
tracking the motion of the sun. The results account for effects 
of solar radiation intensity, exposure time, and incident angle 
of the sun light. Similar to the exposure time, the available 
energy changes from one season to another. Figure 7 shows 
four instances of the surface available energy (fixed array) per 
unit area for one Phobos orbit around Mars at equinoxes and 
solstices. 


Available energy for Tracking Solar Array—represents the 
best scenario for energy generation. It assumed the solar array 
is able to track the motion of the sun perfectly such that the 
normal direction of the solar array is parallel to the sun vector. 
The result accounts for effects of solar radiation intensity and 
exposure time. Figure 8 shows four instances of the surface 
energy (sun tracking array) per unit area for a Phobos orbit at 
equinoxes and solstices. 


4. USAGE EXAMPLE 


This section presents a use case demonstrating the usefulness 
of these results in terms of mission analysis and vehicle 
design. The data is generated using a 450,000 facet STL 


Phobos model. 


Lighting Condition at 11 Science Sites 


Eleven preliminary sites were selected for their scientific 
exploration value. Table 2 lists the longitude and latitude 
coordinate of these sites. 


Figure 9 to Figure 11 present the lighting conditions at these 
11 sites as function of time. 


As expected, exposure time for sites near the equator is 
relatively constant over time compared to the exposure time 
for sites further from the equator. Exposure time for sites 5, 6, 
9, 10, and 11 vary significantly with Martian seasons. Results 
for available energy per unit area for fixed arrays include 
the effects of exposure time, radiation intensity and incident 
angle. Even though exposure time is similar during summer 
and winter for sites near the equator, radiation intensity 1s not. 
The difference in available energy per unit area is significant 
between summer and winter. Compared to the fixed array, a 
sun tracking array could generate more energy per unit area, 
especially for sites in the northern hemisphere such as sites 5, 
6, 9, and 10 during the winter. 


Preliminary Power Subsystem Design Tool 


Making use of the lighting condition results generated for 
the 11 science sites, a preliminary power subsystem analysis 
tool for a Phobos habitation lander was developed. The 
tool can generate battery and solar array sizing information 
based on the power load and solar array efficiency inputs. 
Equations are presented in the Appendix. Input parameters 
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Figure 7. Available energy per unit area for fixed solar array at equinoxes and solstices 
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Figure 8. Available energy per unit area for sun tracking solar array at equinoxes and solstices 
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Table 2. Coordinate of 11 science sites on Phobos 
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Figure 10. Energy per unit area per Phobos’ orbit for a fixed array for site 1 to 11 (kW - hr/m? 
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Figure 11. Energy per unit area per Phobos’ orbit for a sun tracking array for site 1 to 11 (kW - hr/m7?) 


for an example analysis are in Table 3. 


Table 3. Example parameters for a habitation lander 


Attribute 


Solar array total efficiency 20 
Solar array type sun tracking 
Power load 10.0 


Battery max discharge % 80) 





Figure 12 shows the minimum surface area of the solar array 
needed for the lander as a function of time, and Figure 13 
shows the minimum battery capacity needed for the lander as 
function of time, considering the desired maximum discharge 
percentage. 


These results provide preliminary information for engineers 
in designing the power subsystem of the habitation lander 
operating at the 11 science sites. 


5. CONCLUSION 


Many studies have suggested that a human mission to Phobos 
is an important precursor and catalyst for NASA’s human 
missions to Mars [1], [2], [3], [4]. Among other environ- 
mental factors, lighting conditions play a crucial role for 
mission success. This study used a high fidelity computer 
simulation to investigate the lighting conditions on Phobos. 
The simulation includes an environment model, a Phobos 
model, and an occultation model. These models work to- 
gether to provide reliable results of the lighting conditions on 
Phobos. The results show that the lighting condition changes 
significantly from one Martian season to another: radiation 


intensity changes from its minimum during the summer to 
its maximum during the winter; solar eclipse length is long 
during the spring and fall, and no solar eclipse occurs during 
the summer and winter; exposure time per orbit is relatively 
uniform over the surface during the spring and fall but varies 
with latitude during the summer and winter; a sun tracking 
solar array could generate 50 to 75 % more energy than a 
fixed solar array. It is hoped that these results will be useful 
in the vehicle design process for a future mission to Phobos. 


APPENDIX 
Radiation intensity 
1 
I Phobos = Lau pe (1) 
Phobos 


Where J pPpovos 1S the solar radiation intensity at Phobos, J; 4y 
(+1367 W/m2) is the solar radiation intensity of 1 Astro- 
nomical Unit (AU), and dpponos 18 the distance of Phobos 
with respect to the Sun in AU. 


Radiation flux 


& = I Phobos COS (0) Ovapode (2) 


® is the radiation flux for a specific facet. 6 is the angle 
between the Sun vector and a facet’s normal vector. derpose 
is the Kronecker delta function, which equals one when the 
facet is exposed to sun lighting or zero otherwise. 
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Figure 13. Minimum battery capacity needed for habitation lander (kW - hr) 
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Where Agrray 18 the solar array area, tsnadowed 18 the amount 
of time that the facet is in shadow for one Phobos’ orbit, 
Pavg 18 the average power load for the vehicle, 7) is the solar 
array total efficiency, and €gyailable 18 the total available solar 
energy per unit area for one Phobos’ orbit. 


Minimum battery 
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Where Coattery 18 the battery capacity, IMD is maximum 
discharge factor, and the constant coefficient is a correction 
factor to accommodate the additional energy that may be 
needed when power generation is smaller than the power load 
(consuming battery energy). 
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